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Abstract. Using Gabor analysis, we give a complete characterization of all 
lattice sampling and interpolating sequences in the Fock space of polyana- 
lytic functions, displaying a "Nyquist rate" which increases with n, the de- 
gree of polyanaliticity of the space. Such conditions are equivalent to sharp 
lattice density conditions for certain vector- valued Gabor systems, namely su- 
perframes and Gabor super-Riesz sequences with Hermite windows, and in the 
case of superframes they were studied recently by Grochenig and Lyubarskii. 
The proofs of our main results use variations of the Janssen-Ron-Shen duality 
principle and reveal a duality between sampling and interpolation in polyan- 
alytic spaces, and multiple interpolation and sampling in analytic spaces. To 
connect these topics we introduce the polyanalytic Bargmann transform, a 
unitary mapping between vector valued Hilbert spaces and polyanalytic Fock 
spaces, which extends the Bargmann transform to polyanalytic spaces. Mo- 
tivated by this connection, we discuss a vector-valued version of the Gabor 
transform. These ideas have natural applications in the context of multiplex- 
ing of signals. We also point out that a recent result of Balan, Casazza and 
Landau, concerning density of Gabor frames, has important consequences for 
the Grochenig-Lyubarskii conjecture on the density of Gabor frames with Her- 
mite windows. 



1. Introduction 

In this paper we find a new, and perhaps unexpected, connection between polyan- 
alytic functions and time- frequency analysis, and use it to obtain a complete charac- 
terization of all lattice sampling and interpolating sequences in the Bargmann-Fock 
space of polyanalytic functions or, equivalently, of all lattice vector-valued Gabor 
frames and vector- valued Gabor Riesz sequences for L 2 (M.,C n ). 

1.1. Overview. The Bargmann-Fock space of polyanalytic functions, F™(C d ), con- 
sists of all functions satisfying the equation 

(1.D (sfm-o, 

and such that 

(1.2) / \F{z)\ 2 e~^dz < oo. 
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Functions satisfying (II. 1|) are known as polyanalytic functions of order n. Since 
(jl.ip generalizes the Cauchy-Riemann equation 

then the space F n (C d ) is a generalization of the Bargmann-Fock space of analytic 
functions, J-(C d ) = F 1 (C d ). In the case of J-(C), a complete description of the sets 
of sampling and interpolation is known [5H] . |46| . |47] . 

Polyanalytic functions inherit some of the properties of analytic functions, often 
in a nontrivial form. However, as in the theory of several complex variables, many of 
the properties break down once we leave the analytic setting. An obvious difference 
lies in the structure of the zeros. For instance, while nonzero entire functions do not 
have sets of zeros with an accumulation point, polyanalytic functions can vanish 
along closed curves: just take F(z) = zz — 1 = |z| 2 — 1, a polyanalytic function of 
order 2. Polyanalytic functions have been investigated thoroughly, notably by Balk 
and his students [8]. They are naturally related to polyharmonic functions, which 
have an intriguing structure of zero sets [57], [pj, [4"T] . 

We will study the spaces F"(C d ) using time-frequency analysis, and think about 
the polyanalytic functions in a different way from the classical classical approach. 

The link to time-frequency analysis has impressive consequences: by endowing 
the spaces F"(C d ) with the structure inherent to translations and modulations, 
one can use tools that were unavailable with complex variables. This will provide 
F n (C d ) with properties reminiscent of the classical analytic Fock space. 

In order to state our main results, we briefly recall some definitions. See sections 
4 and 5 for more details. We associate the sequence A = {(x,w)} with the sequence 
of complex numbers T — {(x + iw)}. Then T is an interpolating sequence for F n (C d ) 
if, for every sequence {ca,j} & I 2 , there exists F £ F"(C d ) such that 

e* mK "-tWV(«) = o ili , 

for every z £ T. We say that r is a sampling sequence for F n (C d ) if there exist 
A,B >0 such that, for every F G F n (C d ), 

A \\F\\ 2 FHCd) < \F(z)\ 2 e-^' 2 < B ||F|| 2 „ (c .) . 

The concept of interpolating sequences has its roots in deep problems in com- 
plex analysis, and sampling sequences are a major issue in signal processing, since 
they correspond to the case where stable numerical reconstructions of a function 
from its samples is possible. Monograph [35] is a good introduction to sampling 
and interpolation and its interconnections with other branches of pure and applied 
mathematics. 

Our main results, Theorem 4 and Theorem 6 below, use the concept of Beurling 
density, which, in the lattice case, is given by D(T) = D(A) = |detA| - , where 
A = Al?. 

Theorem 4. The lattice T is a sampling sequence for F"(C) if and only if: 

D(T) > n. 

Theorem 6. The lattice T is an interpolating sequence for F n (C) if and only if: 

D(T) < n. 
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These results follow by establishing one duality between sampling in F n (C d ) 
and multiple interpolation in ^"(C 1 ) and a second duality between interpolation in 
F n (C d ) and multiple sampling in ^(C 1 ). When d = 1, these properties allow us to 
directly apply the results in [12J . Taking n — 1 we recover the well known duality 
between sampling and interpolation in .F(C d ). 

Theorems exhibiting a "Nyquist rate" phenomenon tend to be hard to prove. 
They have been studied, for general sequences, in spaces of analytic functions, first 
in the Paley- Wiener space [TT] , [35] , [37] and then in Bargmann-Fock [35] . [46] . [4"T] 
and Bergman [35] spaces of analytic functions. There are two reasons for us to 
believe that Beurling type methods do not work here. First, they use complex 
variables tools that are not available in the polyanalytic situation. Second, in all of 
these spaces, there were known "doubly orthogonal systems" [45], which provided 
the eigenfunctions for the fundamental equation involving the "concentration op- 
erator" . We are not aware of a system with this double orthogonality property in 
the polyanalytic situation. 

Therefore, we introduce new tools. 

With a view to relating our problem to one concerning the density of vector val- 
ued Gabor systems, we extend Bargmann's work [10] to the setting of polyanalytic 
functions. Once we do this, our argument, which is conceptual in nature, follows 
smoothly. 

It is also worth noting that the density Theorem in Gabor analysis has itself a 
very rich story, beginning with fundamental but imprecise statements by John Von 
Neumann and Dennis Gabor, which caught the attention of mathematicians after 
conjectures by Daubechies and Grossman [15]. See the survey article [28], [13] and 
the important special windows studied in [35j . 

To give a context to our approach, recall the connection between the classical 
(analytic) Bargmann-Fock space and time-frequency analysis. 

It is well known that, up to a certain weight, the Gabor transform with a gaussian 
window belongs to the Fock space of analytic functions. Moreover, it has been 
shown that this is the only choice leading to spaces of analytic functions [3J . 

However, a nice picture emerges when we take Hermite functions as windows. 
The analytic situation generated by the gaussian window then becomes the tip of 
the iceberg of a larger structure involving spaces of polyanalytic functions. Indeed, 
the Gabor transform with the nth Hermite function is, up to a certain weight (the 
same as in the analytic case), a polyanalytic function of order n + 1. 

To fully understand the situation, we will need the spaces constituted by the 
functions satisfying (|1.2[) . which are polyanalytic of order n, but are not polyanalytic 
of any lower order (in particular they have no analytic functions). These are the 
true polyanalytic Fock spaces T n (<C d ). The polyanalytic Fock and true polyanalytic 
Fock spaces are related by the following orthogonal decomposition (see Corollary 1 
in section 3): 

F"(c d ) = j=°(c d ) e ... e r 1 - 1 ^). 

Then, each space J 7n (C < *) is associated with Gabor transforms with the nth Her- 
mite window. Such occurrence, which seems to have been hitherto unnoticed, will be 
fundamental in our discussion. This observation is related to some recent develop- 
ments in Gabor analysis with Hermite functions [23] , [24] , [20] , to Janssen's approach 
to the density Theorem [31], [33] and also to the techniques used in [29] . [30] . [53] . 
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which suggest that wavelet spaces and polyanalytic functions share intriguing pat- 
terns. 

Fock spaces of polyanalytic functions are briefly mentioned in Balk's monograph 
[8j and they are implicit in quantum mechanics, in connection with the Landau 
levels of the Schrodinger operator with magnetic field [44] , [21] and displaced Fock 
states [51] . However, we were not able to find any reference to polyanalytic func- 
tions in the mathematical physics literature, apart from [52] . where creation and 
annihilation operators are used. 

1.2. The results of Grochenig-Lyubarskii and of Balan-Casazza-Landau. 

Our results are connected to a very recent result of Grochenig and Lyubarskii, 
which deserves more specific comment in this introduction. Denote by G(h„,A) 
the set of the translations and modulations indexed by the lattice A and acting 
coordinate-wise on the vector h„. Then, the following result holds. 

Theorem [24] : Let h„ = (ho, h n — \) be the vector of the first n Hermite 
functions. Then G(h n ,A) is a frame for L 2 (M., C") if and only if 

D(A) > n. 

One may wonder if the equivalence of this condition to the one in Theorem 4 reflects 
causality or casuality. As we will see, causality is the answer. Actually, one of the 
key steps in our approach consists of showing that Theorem 4 is equivalent to the 
above Theorem. 

The original proof of this Theorem in [21] combines the use of the so-called 
Wexler-Raz biorthogonality relations with complex analysis techniques based on 
properties of the Weierstrass sigma function. We will give an alternative proof, 
which is considerably shorter (at the cost of using deeper results from the litera- 
ture) and has the advantage of also characterizing the vector valued Gabor Riesz 
sequences with Hermite windows in the following statement, which is equivalent to 
Theorem 6: 

Theorem 7: Q(h n ,A) is a Riesz sequence for L 2 (M, C") if and only if 

D(T) < n. 

One should also remark that the necessity of the condition D(A) > n for vector 
valued frames, follows from a result of Balan [5]. Moreover, this condition holds for 
general sets, since a Ramanathan-Steger [40] type argument is used. 

Now let us look closer at the problem of deciding when the translations and 
modulations of a single Hermite function constitute a frame. 

It can be easily seen as a corollary of the above Theorem that D(A) > n is suffi- 
cient for the system G(/i„_i, A) to be a frame. It has been observed by Grochenig 
and Lyubarskii that there are some examples which support the intriguing conjec- 
ture that such a result might be sharp. 

A recent result of Balan, Casazza and Landau shows that this conjecture cannot 
hold for general sets. Let So stand for the Feichtinger algebra [16] . 

Theorem [7] Assume G(g,A) is a Gabor frame for L 2 (R d ) with g £ Sq. Then, 
for every e > ; there exists a subset J e C A so that G(g, J e ) is a Gabor frame for 
L 2 (R d ) and its upper Beurling density satisfies D + (J e ) < 1 + e. 
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This Theorem implies that Grochcnig-Lyubarskii conjecture can only be true for 
lattices. Indeed, for every n, the So norm of the Hermite functions (see [34]) is 



"TV. 

and therefore, 

h n e So- 

Thus, the Balan-Casazza-Landau Theorem tells us that we can always find, for every 
n, a set J with Beurling density arbitrarily close to one and such that G(h n ,A) is 
a frame. However, it is still possible that the conjecture is true for lattices, since a 
subset of a lattice may not be a lattice. 

Conjecture 1. Let A be a lattice mffi 2 . If G(h n ,A) is s Gabor frame for L 2 (R d ), 
then 

D{A) > n. 

1.3. New concepts. We will introduce some new concepts troughout this paper, 
in order to extend Bargmann's theory [10] and its connection to Gabor analysis in 
the polyanalytic setting. 

• We first introduce what we call the true-poly analytic Bargmann transform: 

(B n f)(z) = (TrWnO-ie'W [e^ F(z) 

Here F stands for the Bargmann transform of /. As we will see, this is a 
unitary mapping from L 2 (K. d ) to !F n (C d ). This mapping relates to Gabor 
transforms with Hermite windows $„ in the following way: 

V*J(x,u) = e i ** u -* ]s £-{B n f){z), 

and we will provide its basic theory starting from this relation, following 
as much as possible the presentation of the Bargmann transform given in 
Grochenig's book [HI section 3.3]. 

• For vector- valued functions f = (/o, /n-i), we define the polyanalytic 
Bargmann transform, 

(B"f)= y, ( Bfe /*)> 

0<fc<M-l 

which will be unitary between L 2 (R d , C") and F n (C d ). 

• In the last section we will see that the polyanalytic Bargmann transform is 
a special case of a vector- valued version of the Gabor transform. Although 
this transform plays no direct role in the proofs of the main Theorems, 
it must be in the picture for completeness, since it is the natural time- 
frequency transformation of which the polyanalytic Bargmann transform is 
a special case: 

n-l 

V g i(x,uj) = J2 V 9kfk(x,uj). 

k=0 

In the case where {<?fc}fc=o,...n-i constitutes an orthonormal sequence, it 
provides an isometry between L 2 (R d , C n ) and between L 2 (R 2d ). This trans- 
form is the continuous counterpart of the superframes discussed in [5B] , [3] , 
[24] Theorem 2.7]. 
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1.4. Technical summary of proofs. With the tools described above at hand, 
our main argument will depend on two profound results. More specifically, we 
will combine variations on the Janssen-Ron-Shen duality principle [43j with the 
characterization of multiple sampling and interpolation sequences in the Fock space 
[12] . The duality principles reflect all the rich inner structure of Gabor frames. 
The second result uses a deep elaboration on Beur ling's balayage technique [TTj 
developed by Seip in [48] . 

We will proceed as follows. 

First, using an orthogonal basis for the polyanalytic Fock spaces, we prove the 
unitarity of B n and B™. Then we study sampling in F"(C). Using the unitary 
mapping B™ we show that the problem is equivalent to the study of vector valued 
frames with Hermite windows, also known as superframes [4] . [24] . This problem 
has been recently studied in [24] , but we provide an alternative proof, which is more 
natural in the context of sampling and interpolation: applying a vector valued ver- 
sion of the Janssen-Ron-Shen duality we translate the statement into a problem 
concerning unions of Riesz sequences. After noticing that the latter is equivalent 
to a multiple interpolation problem in Fock spaces of analytic functions, we apply 
the interpolation result in [12]. We then study interpolation in F"(C). In order 
to do this, we "dualize" the arguments that we have used in the sampling part, 
once again using the Janssen-Ron-Shen duality, this time between vector-valued 
Riesz sequences and multi-frames with Hermite functions. This translates our in- 
terpolation problem into one of multiple sampling. Noticing that this problem is 
equivalent to multiple sampling in Fock spaces, we apply the sampling result from 

HZ|. 

1.5. Organization of the paper. The next section contains the classical tools 
that we are going to use. We list the basic properties of the Gabor transform, the 
Bargmann transform and the Hermite functions. 

In the third section, we introduce the true polyanalytic Bargmann and the poly- 
analytic Bargmann transforms. By making a connection with the Gabor transform, 
we study their basic properties, find an orthogonal basis for the polyanalytic Fock 
spaces and prove the unitarity properties. 

Our main results are in the fourth and fifth sections, where we derive the duality 
principles and study sampling and interpolation for F™(C). 

In our last section we introduce the super Gabor transform and make some 
remarks of a more informal character concerning applications and open problems. 

2. Background 

2.1. The Gabor transform. Fix a function g ^ 0. Then the Gabor (short-time) 
Fourier transform of a function / with respect to the "window" g is defined, for 
every x, u> € M. d , as 

(2.1) V g f(x,u)= [ f{t)W^)e- 2 ^dt. 

The following relations arc usually called the orthogonal relations for the short-time 
Fourier transform. Let /i,/ 2 ,ffi,ff 2 € L 2 (R d ). Then V g J ll V g2 f 2 € L 2 {R 2d ) and 



(2.2) 



(Vg 1 fl,Vg 2 f2) L 2( M 2d) - (fl,fc) L *(R*) (gi,92) L 2(^dy 
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The Gabor transform provides an isometry 

V g : L 2 (R d ) -> L 2 (R 2d ), 

that is, if f,g G i 2 (M d ), then 

( 2 -3) ||Vg/|| i2 ( R 2d) = ||/Hl 2 (r< 1 ) ll3llL 2 (R d ) ■ 

For every x, u> G K d define the operators translation by a; and modulation by ui as 

r s /(t) = f(t-x), 

M u f{t) = e 2 ^ */(*)• 
Using these operators we can write (|2.1j) as 

y g /(x,cj) = (f,M u T x g) L2(ud) . 

2.2. The Bargmann transform. Here we will use multi-index notation: z = 
(zi,...Zd), n = (m,...,nd) and |n| = ni + ... + n^. The Bargmann transform, 
defined by 

(Bf)(z)= [ f{t)e Mz ^-^dt, 



is an isomorphism 

B:L 2 {R d )^T{<£ d ), 

where F(C d ) stands for the Bargmann-Fock space of analytic functions in C d with 
the norm 

(2.4) \\F\\% m = / \F(z)\ 2 e-^dz. 

Jc d 

The collection of the monomials of the form 

(2-5) e n (z)= — z» = JJ -=z% 

where n — (rii^.nd), with rii > 0, constitutes an orthonormal basis of !F(C d ). 
The reproducing kernel of J-(C d ) is the function e™ '. This means that, for every 
F G JF(C d ), 

(F( II ),e'»}, (ci) =F(4 
Differentiating n — k times the corresponding reproducing equation, we obtain 
(2.6) (F(w),w n - k e^ w ) ncd) = n k - n F^ k \z), 

where for d > 1 we are using the multi-index derivative 



c?z dz\...dzd 

A simple calculation shows that the Bargmann transform is related to the Gabor 
transform with the Gaussian window tp(t) = 2ie~ 7rt by the formula 

(2.7) V v f(x, -ui) - e " ^ (Bf)(z), 
where z = x + iu. 

We will need one more operator. Define a "translation" fit on !F(C d ) by 

(2.8) p z F(C) = e^-^e^Cp^ _ z y 
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The operator j3 z satisfies the intertwining property 

(2.9) P Z B = BM U T X , z = x + iuj. 

2.3. The Hermite functions. The Hermite functions can be defined via the so 
called Rodrigues Formula 

MO-**- 1 (f )>"*)• 

where c n is chosen in such a way that they can provide an orthonormal basis of 
i 2 (M). Now let n = (n\, ■■■ : rid) and x S M d . The d- dimensional Hermite functions 
are 

*n(a>) = Il^ife)' 

3=1 

They form a complete orthonormal system of L 2 (IR d ). 

A very important property of the Hermite functions (see for instance |34j ) is 
that they are mapped onto a basis of the Bargmann-Fock space via the Bargmann 
transform. 

(2.10) (B*n)(«) = e n (z). 



3. POLYANALYTIC FOCK SPACES AND POLYANALYTIC BARGMANN TRANSFORMS 

3.1. Definitions. In this section we use multi- index notation in such a way that 
there will be no difference between the one and the d-dimensional case. Only at 
the end of the last two sections is it necessary to specialize d — 1 . 

It is well known [5] that every polyanalytic function of order n can be uniquely 
expressed in the form 

(3-1) F(z)= Yl *<Pk(*), 

0<k<n-l 

where {<£>p(z)}p=o are analytic functions, each of them with a power series expan- 
sion, 

j>0 

As a result, there is also a power series expansion for the polyanalytic function F: 
(3.2) F(z)= ^' 

0<p<n-l j>0 

We will often use the inner product in the polyanalytic Fock space, given by 

(F,G) Fn(cd) = [ F{z)G{z)^ 2 dz. 
Observe also that this implies 

(F,G) Fn , cd) = (e-^F,e-^G) 

\ I L 2 (R 2d ) 
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3.2. The true polyanalytic Bargmann transform. 

Definition 1. The true polyanalytic Bargmann transform of order n, of a function 
on R d , is defined by the formula 

rf n r 

(3.3) {B n f)(z) = frM n \)-*e*W — [e~^ F(i 
where F{z) = (Bf)(z). 

Clearly B°f = Bf and B n is a generalization of the Bargmann transform. We 
now provide the fundamental properties of B n . We try to stay as close as possible 
to the presentation of section 3.4 in [22 . The next proposition is the departing 
point of our study. 

Proposition 1. If f is a function on M. d with polynomial growth, then its true 
polyanalytic Bargmann transform B n f is a polyanalytic function of order n + 1 on 
C d . If we write z = x + iui, then this transform is related to the Gabor transform 
with Hermite windows in the following way: 

(3.4) V 9n f(x,u) = e l ' KXUJ - 1T ^(B n f)(z). 
Moreover, if f £ L 2 (R), then 

(3-5) Ip /Ilz2(c<*,e-"M 2 ) = ||/|| L 2( R d) • 

Proof. Let F — Bf. The following calculation is from Proposition 3.2 in [53], where 
(HHD is used: 

V# n f(x,w) = (/, A4T^$„) L2(R£i) 

= {Ff3 z B^ n ) Hcd) 

= (Tr^nirh™™-^ 2 (F(w),e™ w (w-zr) Hcd) 

= (Tr^n!)-^"™^-^ 2 ! 2 ^ Cf\(-7rz) k F^- k \z). 

0<k<n ^ ' 

Now, since the Bargmann transform of a function in M. d is an entire function [2H 
Proposition 3.4.1], the functions F^ n ~ k \z) are also entire, and from (|3.ip we rec- 
ognize the sum as a polyanalytic function of order n + 1 . To prove (|3.4|) observe 
that the last expression can be written as 



.\A* 



giTrxw-TT-a-^lnl^-igjr^l " F{z)\ = j™"-*^ (B n f)(z). 

The isometric property p.5p is an immediate consequence of (|3.4[) and (|2.3p . □ 
3.3. Orthogonal decomposition. 

Definition 2. For k,m G Nq, consider the functions ek, m defined as 

(3.6) e k , m (z) = (Tr^kiyh^ 2 [e~^ e m {z) . 

From (|2.10[) one can easily see that 

(3.7) e k , m {z) = (£ fc $ m )(z). 

Proposition 2. The set {efc,m} 0<fe< „_ 1 . m >o* s an orthogonal basis of F n (C d ). 
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Proof. The orthogonality follows from (l3J|) . (13T4j) and (j2T2[> . since 



To prove completeness of {e^m} in F"(C d ), suppose that F G F ra (C d ) is such that 
(F, ek, m )yr (C d) = 0, 0<fc<n-l; to > 0. 



For fc = 0, we can use the representation of F in power series (I3.2[) . Interchanging 
the sums with the integrals and using the orthogonality of the functions (12. 5p , the 
result is 

(3.8) (F,e , m ) ncd) = <^™'V /^LZlrn\ = °' m -°' 

For fc > 1, a calculation using integration by parts gives: 



{F, e k , m ) jr (cd) = / F(z)e k , m (z)e *< z ' dz 

= I e-^ 2 7^)p...(p-k + l) V zP^Y^c^dz 

J ^ k<p<n-l j>0 



k<p<n-l j>0 



As a result, 



Ep...(p — k + l)(p + m — k)\ 
1 —. ^ % i= L c m+P -k, P = 0, m > 0, < k < n - 1, 

fc<p<n— 1 

resulting in a triangular system for each to. Solving this system we obtain Cj tP = 
for k < p < n — 1 and j > 0. Therefore, F = 0. □ 

Remark 1. Clearly the orthogonality in Proposition 2 can be obtained directly by 
integration by parts and moving to polar coordinates. For k = this has the ad- 
vantage of showing that the functions are also orthogonal in the polydisk, providing 
the useful "double orthogonality" property as in the proof of |22[ Theorem 3.4.2.]. 
However, for k > 0, the boundary behavior required in the integration by parts 
eliminates this advantage, making the functions ek >m less likely to possess such a 
property. 

Remark 2. It is clear that these functions are reminiscent of the so-called special 
Hermite functions, which are the Wigner transforms of two Hermite functions [50j . 
They also appear in the study of Landau levels in [21j . 

Definition 3. The true poly analytic Fock space of order n is defined as 
(3.9) r l (C d ) = Span [{e n , m (z)} m>a 
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Remark 3. Observe that 



rl \ " r i 1 rfm+n 



dz 



dz n dz n 



e 



Therefore, our functions e n , m are essentially the complex Hermitian functions in- 
troduced in |44l pag. 126] and, as a result, according to Theorem 7.1 in [44], the 
true poly analytic Fock spaces are the eigenspaces of the Schrddinger operator with 
magnetic field in M 2 , associated with the eigenvalue n + |. Also, observe that the 
basis used in |39j approaches this one by a formal limit procedure. 

The orthogonal basis property has the following consequence: 

Corollary 1. The polyanalytic Fock space, F n (C d ), admits the following decompo- 
sition in terms of true polyanalytic Fock spaces T (C ). 

(3.10) F"(C d ) = .F°(C d )©...©.F n - 1 (C d ). 

This results in a definition equivalent to the one in [53j . where the spaces were 
defined using the decomposition. Observe that F 1 ^) = T°{C d ) = T(C d ) and 
that functions in T n (<C d ) are polyanalytic of order n + 1. 

3.4. Unitarity of B n . The true polyanalytic Bargmann transform keeps the uni- 
tarity property 

Theorem 1. The true polyanalytic Bargmann transform is an isometric isomor- 
phism 

B n : L 2 (R d ) — > T n {€- d ). 



Proof. Since we know from (|3.5|) that B n is isometric, we only need to show that 
B n [L 2 (M. d )] is dense in JF™(C rf ). This is now easy, since the Hcrmitc functions 
constitute a basis of L 2 (R d ) and, by (|3.7|) . they are mapped onto the basis {e n , m (z)} 
of «F"(C). Since B n [L 2 {R d )) contains a basis of F l (C d ), then it must be dense. □ 

3.5. The polyanalytic Bargmann transform. Now consider the Hilbert space 
H = L 2 (R d ,C n ) consisting of vector-valued functions f = (/o> /n-i) with the 
inner product 



(3-11) (f,g) W = (fk,9k) L 2 {Rd y 

0<fe<n-l 

The polyanalytic Bargmann transform of a function f = (/o, f n -i) is defined 

as 

(3.12) (B"f)(z)= E (B k f k )(z). 

0<fc<n-l 

The next Theorem, which may have independent interest as a generalization of 
Bargmann 's unitary transform, will be the cornerstone in the proof of our main 
results regarding sampling and interpolation. 

Theorem 2. The polyanalytic Bargmann transform is an isometric isomorphism 

B" : H -> F n (C d ). 
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Proof. For the isometry, first observe that, from (|2.2p and (|3.4p . 

(B k fk,B' J fj) rn{cd) = S k ,j. 
Then, using the isometric property of B n , 



||B n f||p„ (C£i) - ^2 ||^ fc /fe||^ (C d) 

0<k<n-l 



^ ll/fellL 2 ^) — \\f\\ 



2 

\n ■ 

0<k<n-l 



Moreover, B n [L 2 (R d )} is dense in F"(C d ), since, by the decomposition (|3.10p . every 
element F € F™(C d ) can be written as 

F= £ ^ 

0<fc<n-l 

with F fe S T k (C d ), < ft < n - 1. Since Z? fe is unitary, there exists / fc G J L 2 (K d ) 
such that F/j = B k ft, for every < ft < n — 1. It follows that F = B n f, with 

f =(/ ,...,/n-l). □ 

4. Sampling in F"(C) 

4.1. Definitions. We will work with lattices. A lattice is a discrete subgroup in 
R 2d of the form A = AZ 2d , where A is an invertible 2d x 2c? matrix. We will define 
the density of the lattice by 

D(A) - W 

If we write z — x + iu> and 7r z g = M u T x g, the adjoint lattice A is defined by the 
commuting property as 

A = {(it <E R 2d : 71-271"^ = tt^tTz, for all z e A}. 

If A = aZ x /3Z, then A = ^ _1 Z x a -1 Z. In general, 

A = J{A T Y lr L d , 



I 

-I 



(consisting of d x d blocks) is 



where A T is the transpose of A and J = 
the matrix defining the standard sympletic form (see [24] and [17]). Therefore, 

(4-2) D(A°) = -4tt- 

17 v 7 D(A) 

We will use the notation T = {z = x + iuj} to indicate the complex sequence 
associated with the sequence A = (x,oj). The density of T will be the density of 
the associated lattice, that is D(T) = D(A). 

Definition 4. T is a sampling sequence for the space F™(C d ) if there exist positive 
constants A and B such that, for every F 6 F"(C rf ), 

(4-3) A \\Ff Fn{cd) < J2 \F(z)\ 2 e~^ 2 < B \\F\\ 2 Fn(cd) . 

zer 
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The definition of sampling in the spaces !F k (C ) is exactly the same. 

Now we take the following definition, obtained from [12j page 114], by making 
a small simplification (in the notation of [T2] page 114] we set v(z) — n ) and 
rewriting it in our context (observe that the weight e l7TXUJ makes no diference). 

Definition 5. A sequence T n , consisting of n copies ofT is a multiple interpolating 
sequence in the Fock space J-(<C d ) if, for every sequence {ct^j }fc=o,...n-i such that 

{ a \ k j}k=o,...n-i G I 2 , there exists F G J-(C d ) such that (F,(3 z ek) = c*t] , f° r a H 
< k < n — 1 and every zeT. 

Consider again the Hilbert space Ti = L 2 (M. d ,C n ) consisting of vector- valued 
functions f = (/o, •••! /n-i) with the inner product p. lip . The time-frequency 
shifts act coordinate-wise in an obvious way. 

Definition 6. The vector valued system C/(g,A) = {Mo;Txg}(x,w)eA * s a Gabor 
superframe for Ti if there exist constants A and B such that, for every f G Ti, 

(4.4) M\n 2 H < E i(M« g ) w i 2 <b\u 2 h . 

(i,io)EA 

Superframes were introduced in a more abstract form in [26) and in the context 
of " multiplexing" in [1]. We will need a fundamental structure Theorem from time- 
frequency analysis, namely the following version of the Janssen-Ron-Shen duality 
principle [H Theorem 2.7]. 

Theorem A. Let g = (go, g n ~i)- The vector valued system G(g, A) is a Gabor 
superframe for Ti if and only if the union of Gabor systems U^~Q^(.gfe, A ) is a Riesz 
sequence for L 2 (R d ). 

4.2. Duality principle. In this section we will obtain the following duality prin- 
ciple. 

Theorem 3. T is a sampling sequence for F" (C d ) if and only if the adjoint sequence 
is a multiple interpolating sequence in the Fock space T(jC d ). 

We first prove two Lemmas. Combining them with Theorem A gives Theorem 

3. 

Lemma 1. The union of Gabor systems ^ZoS(hk, A) is a Riesz sequence for 
L 2 (M. d ) if and only if T n is a multiple interpolating sequence in the Fock space 
T(C d ). 

Proof. The union of Gabor systems U^q^/i/c, A) is a Riesz sequence for L 2 (R d ) 

if, for every sequence {a[ y }k=o,...n-i G I 2 , there exists a / G L 2 (R d ) such that 

(/) M u T x hk) — ol\ j, for all k = 0, ...n — 1 and every (x, lu) G A. Using the unitarity 
of B and the intertwining property (|2.9[) gives 

{f,M u T x h k ) = (Bf,l3 z e k ), 

and setting F = Bf shows that r„ is a multiple interpolating sequence in the Fock 
space T{C d ). □ 

The next Lemma is a key step in our argument and it is at this point that the 
unitarity of the polyanalytic Bargmann transform is essential. 
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Lemma 2. Let h„ = (ho, ...,h n -i). Then the set Q(h. n ,A) is a Gabor superframe 
for H. = L 2 (M. d , C n ) if and only if the associated complex sequence T is a sampling 
sequence for F n (C d ). 

Proof. Using the definition of the inner product (|3. 1 1|) , identity (|3.4|) and the defi- 
nition of the polyanalytic Bargmann transform, it is clear that 

(4.5) (f,M w T x h n ) n = Y, (fk,M„T x h k ) LHMd) 

0<fc<n-l 

Y e inx "-%M 2 (B k fk)(z) 

0<fc<n-l 

(4.6) = e i7ra;w ~^ |z|2 (B n f)(z). 

Therefore, setting F — B n f, the unitarity of B" shows that the inequalities (|4.4|) 
are equivalent to (|4.3|) . □ 

4.3. Main result. We will need the concept of Beurling density of a sequence. 

Let I be a compact set of measure 1 in the complex plane and let n~(r) denote 
the smallest (and n + (r) the biggest) number of points from T to be found in a 
translate of rl. We define the lower and the upper Beurling density of T to be 

n (r^\ n~^~(r} 
D~(T) = lim sup — ^ and D + (T) = lim sup — 

r — >oo t T— *oo T 

respectively. When V is associated with the lattice A, D~(T) = D+(T) = D(T) = 
D(A). 

We will now use the following result, which is Theorem 2.2 in [12] . Observe that 
we can remove the uniformly discrete condition from the statement in |12j since we 
are dealing only with lattices. 



Theorem B. The sequence T n is a multiple interpolating lattice sequence in the 
Fock space .F(C) if and only if D(T n ) < 1. 



From this we obtain the characterization of sampling lattices in F™(C). 

Theorem 4. The lattice T is a sampling sequence for F"(C) if and only if D(T) > 
n. 

Proof. We know by the duality principle that T is a sampling sequence for F"(C) 
if and only if the adjoint sequence T° is a multiple interpolating sequence in the 
Fock space .F(C) . By definition of Beurling density, it is obvious that 

D(T° n )=nD(r ). 

Therefore, Theorem B states that T° is a multiple interpolating sequence in the 
Fock space !F(C) if and only if D(T°) < -. Using (|4.2p . we conclude that T is a 
sampling sequence for F"(C) if and only if D(T) > n. □ 

Using Lemma 1, we recover Theorem 1.1 of [24] . 

Corollary 2. Let h„ ~ (ho, h n -\). Then the set Q(h n ,A) is a Gabor super 
frame for H = L 2 (R, C n ) if and only if D(T) > n. 
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5. INTERPOLATION IN F n (C) 

5.1. Definitions. 

Definition 7. The sequence T is an interpolating sequence for F n (C d ) if, for every 
sequence E I 2 , there exists F G F n (C d ) such that e l7rxw -f W F(z) = Ojj, /or 

every zeT. 

Definition 8. T/ie sequence T n , consisting of n copies of T is is said to be a 
multiple sampling sequence for J-(C d ) if there exist numbers A and B such that 

(5.1) A\\F\\ 2 ncd) <Y, E m(3 z e k )\ 2 <B\\F\\ 2 Hcd) . 

zero<fc<n-i 

Definition 9. The set Lf k l ~ Q(gk, A) is said to generate a Gabor multi-frame in 
L 2 (R d ) if there exist constants A and B such that, for every f E L 2 (M. d ), 

(5.2) A||/||^ (ffid) < £ \(f> M « T *9k) L * m \ 2 <B\\f\\ 2 L2{Rd) . 

(i,u)£A0<Kn-l 

The dual of the duality principle contained in Theorem A is now required. As 
stated at the end of [25], it reads as follows: 



Theorem C. The set C?(g, A) is a Riesz sequence for L 2 (R d ) if and only if U^Q(g k , A ) 
is a Gabor multi-frame in L 2 (R d ). 

5.2. Duality principle. Now we prove the following duality. 

Theorem 5. The sequence T is an interpolating sequence for F n (C d ) if and only 
z/r" is a multiple sampling sequence for J-(C d ). 

As in the sampling section, we first prove two Lemmas which, combined with 
Theorem C, give the result. The next Lemma requires only the unitarity of the 
Bargmann transform. 

Lemma 3. The set U^^/ifc, A) is a Gabor multi-frame in L 2 (R ) if and only if 
r„ is a multiple sampling sequence for J-(C d ). 

Proof. Similar to Lemma 1: using the unitarity of B and the intertwining property 
(|2.9p gives (f,M u T x hk) = (Bf,f3 z e k ); setting F = Bf it follows from the unitarity 
of the Bargmann transform that (|5.1[) and (|5.2p are equivalent. □ 

Again, we make the key connection in the next step, where the unitarity of the 
polyanalytic Bargmann transform is required. 

Lemma 4. The sequence T is an interpolating sequence for F™(C d ) if and only if 
G(h ni A) is a Riesz sequence for Ti. 

Proof. The sequence T is an interpolating sequence for F n (C d ) if, for every se- 
quence {a itj } € I 2 , there exists F £ F n (C d ) such that e «™^f \ z \ 2 F{z) = a i;j) for 
every z E T. Using the unitarity of B", we find, for every F £ F™(C d ), a vector 
valued function f e H such that B"f = F or, by (|43 ]l -(|4T6 ]l . (f , M w T x h n ) n = 
giTrxw-f \z\ Therefore, the first assertion is equivalent to say that, for every se- 
quence {cti.j} € I 2 , there exists & f £ H such that (f , M u T x h. n ) H = a.i_j, for every 
z E T. This says that £/(h„, A) is a Riesz sequence for H. □ 
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5.3. Main result. Wc will need the following result, which is contained in Theorem 
2.1 in pj]: 

Theorem D. The sequence T n is a multiple interpolating sequence in the Fock 
space T{C) if and only if D(T n ) < 1. 



As before, we can obtain our main result from this one. 
Theorem 6. The lattice T is an interpolating sequence for F n (C) if and only if 



Proof. We know from the duality principle that T is an interpolating sequence for 
F ra (C) if and only if T n is a multiple sampling sequence for !F(C). Once again 
we have -D(T") = nD(T°). Therefore, Theorem D states that T° is a multiple 
interpolating sequence in the Fock space .F(C) if and only if D(T°) > i. As in 
Theorem 5 it follows that T is an interpolating sequence for F"(C) if and only if 



From this and Lemma 4 we obtain a new result characterizing all the lattices 
which generate vector valued Gabor Riesz sequences with Hermite functions. This 
reveals, at least for lattices, the existence of a critical density for vector-valued 
Gabor systems with Hermite functions. 

Theorem 7. Q(h. n , A) is a Riesz sequence for H. if and only if D(T) < n. 

Remark 4. We should remark that the reason we did not care about the Bessel 
condition in the equivalence of the Riesz sequence and interpolating property, used 
several times in the previous section, is that the Hermite functions belong to Fe- 
ichtinger 's algebra Sq (see [T5].|16j ): 



where ip is the L 2 -normalized Gaussian. As a result they satisfy the Bessel condition 



6.1. The super Gabor transform. The polyanalytic Bargmann transform is an 
instance of a more general transform. Although it plays no role in the proofs of our 
main results, we briefly describe it here for completeness of the picture. 

The super Gabor transform of a function f with respect to the "window" g = (go, ...g n -i) 
is defined, for every x,oj S M. d , as 



D(T) < n. 



D(T) < n. 



□ 




[281 theorem 12]. 



6. Generalizations, applications and open problems 



(6.1) 



VgfCar.w) = (f,M u T x g) 




0<fc<n-l 



That is to say, 



0<fc<n-l 



This defines a map 



V g f :H -^L 2 (R 2d ). 
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In the case when the vector g is extracted from an orthogonal sequence {<?fc}fc>o, 
the essential properties of the Gabor transform are kept. As an example of how 
the results concerning Gabor analysis can be generalized to this setting, we obtain 
the isometric properties and the orthogonality relations (the latter valid under the 
slightly weaker condition of biorthogonality) . 

Proposition 3. If 

( 6 - 2 ) (9i,9j) L 2Qi*)= s i,3> 
then V g f is an isometry between Hilbert spaces, that is 



\Vgf || L 2( R 2d) — ||f ■ 



Proof. Using ([6T2]) and (f2T3]> gives 



|V ff f||^ 2(R2d) = ^ ( V g k fk,V g Jk) L 2 { 

0<fc<n-l 

||Vg |c /fc|| L 2( R 2 £ i) 

0<fc<n-l 

= H/ fe lli 2 (R d ) 
0<fc<?i-l 

= llfllw- 



□ 



The orthogonality relations are valid under the slightly weaker condition of 
biorthogonality. 

Proposition 4. If 

( 6 - 3 ) (ffi,i,5 , 2,i) L 2 (Rti ) = Si,j, 

then V g f satisfies 

(6-4) {V gi J 1 ,V a2 J 2 ) L2{R2d) ^(i 1 ,{ 2 ) n . 
Proof. Using and (TO]) gives 

( V ffl,i f l. V 92, J f 2) L 2 (R 2d ) = X! ( f l.fe 1 f 2,fc) L 2 (R 2d) = (fi,f 2 ) 



H ■ 

0<k<n-l 



□ 



Clearly, when we take g = (<&o, ■■■®m-i), we have the following relation with the 
polyanalytic Bargmann transform: 

(B"f)(z) = e inxu - wl ^y s f(x,u;). 

This observation removes some of the mystery from the previous sections. Now we 
are in a position to say that the polyanalytic Bargmann-Fock spaces play exactly 
the same role as the Bargmann-Fock space in the scalar case. We have thus all the 
basic ingredients to build a theory of vector valued (or super) Gabor analysis: 

• A vector valued Gabor transform. 

• A discrete theory for L 2 (R d ) frames and Riesz basis. 

• A special vector of windows providing the connection with complex analysis, 
where, in the case d = 1, a Nyquist rate phenomenon can be observed. 
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The analyzing vector can be extracted from orthogonal systems other than the 
Hermite functions, though they probably do not lead to very structured situations. 
As a first example we may think of the Haar basis. Other wavelets may also be 
used, but we will not pursue this question further in this paper. 

6.2. Applications. Although we are here dealing mostly with questions of a con- 
ceptual nature, there are potential applications of these results in multiplexing, an 
important method in telecommunications, computer networks and digital video, as 
indicated in 0] and )24| . The idea of multiplexing is to encode n independent signals 
fk S L 2 (R d ) as a single sequence f that captures the time- frequency information 
of each one. With suitable windows g = (go, ...g„_i), the time-frequency content 
of each signal can be measured by the associated super-Gabor systems. The super 
Gabor transform V g f(x, w), gives the precise value one wants to approximate via 
the discrete systems. 

6.3. Further questions. As in the scalar-valued case, when the connection to 
complex analysis is missing the lattices generating super-frames should be hard to 
describe as, for instance, the case of a "rectangular" window. What we mean by 
rectangular window is one obtained from the Haar basis, which is the natural gen- 
eralization of the characteristic function of an interval. We wonder how Janssen's 
tie [32] would generalize to this situation. 

It is still unclear if there is a relation between the Bargmann-Fock space of 
polyanalytic functions and vector valued coherent states, as considered in [1] and 
0- 

A rather mysterious topic is sampling and interpolation in the true polyanalytic 
space J-"" _1 (C). Partial results follow from ours. Using decomposition (|3.10p . it is 
easy to see that we obtain two propositions concerning J 7 " (C): one, a sufficient 
condition for sampling, the other a necessary condition for interpolation. 

Proposition 5. If D(A) > n then T is a sampling sequence for .T 7 ™ (C). 

Proof. This is obvious from decomposition ()3.10|) because, if D(A) > n, then in- 
equality (|4.3p holds for every F £ JF"(C). In particular it also holds for every 

FeF"(C). □ 

Proposition 6. IfT is an interpolating sequence for J 7n ^ 1 (C), then 
(6.5) D(A) < n. 

Proof. If D(A) > n, then T is not an interpolating sequence for F" (C). As a result, 
there exists {cti.j} £ I 2 , such that it is impossible to find F £ F n (C), verifying 
giTi-xw-f \z\ p( z ^ — oti,j- Again, from the decomposition (|3.10p . one sees that in 
particular it is impossible to find F £ JF"(C) verifying e mxu -i\A F(z) — ctij. □ 

We may wonder whether these conditions are sharp. 

In the case of Proposition 6 the answer is a definite " no" . This is because, if 
r is an interpolating sequence for ,F n (C), then {(/, M^Txhn-i) L2 ( R )}( x .w)eA is 
a Riesz basis for L 2 (M). As a result, {(/, MuTxhn-i) L2 m\}(x,w)eA° is a frame for 
L 2 (M.) and consequently, using Rieffel-Ramanathan-Steger Theorem [42], 00], we 
must have D(A°) > 1 and, by scalar Ron-Shen duality, this implies D(A) < 1. 
Therefore, estimate (|6.5j) is far from being sharp. 

Let us look closer at Proposition 5. Here things become rather intriguing. 
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It is easy to see that the set Q(h n ,A) is a Gabor frame if and only if the asso- 
ciated sequence T is a sampling sequence for T n {C): since (/, Mu,T x h n -i) L 2, K \ — 
Vh n f(x,w), then (|3.4|) can be written as 

(6.6) (f,MM n ^) L2{R) = e^-^ 2 B n f(z). 

By setting F = B n f, the isometry of B n : L 2 (R) -> F n (C) and the relation (JSH) 
show that definition (|4.3|) is equivalent to the fact that {(/, M^T x h n -i) L2 ^}( x ,w)eA 
is a frame. 

Therefore, Proposition 5 is equivalent to the sufficient condition obtained in [23j . 
where the authors prove that, if D(A) > n, then {(f,M UJ T x h n } L 2^}(x.w)GA is a 
frame and give some evidence to support their conjecture that the result is sharp. 
If true, this would be quite a surprising statement, in face of decomposition (|3.10[) . 
Moreover, by duality, it would bring estimate (|6 . 5[) down to D(A) < l/n. 
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